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Abstract 

We construct a generalized multiplicative renormalization group transforma- 
tion to study the low energy dynamics of a heavy particle tunneling among M 
different positions and interacting with Nf independent conduction electron 
channels. Using a l/A'^y-expansion we show that this M-level scales towards a 
fixed point equivalent to the Nj channel SU{M) x SU{Nf) Coqblin-Schrieffer 
model. Solving numerically the scaling equations we find that a realistic M- 
level system scales close to this fixed point (FP) and its Kondo temperature 

is in the experimentally observable range 1 — WK. 
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One of the simplest examples of non- Fermi-liquid (NFL) systems is provided by the two- 
channel Kondo model where the two conduction electron channels overcompensate the 
spin (or some quasi-spin) of an impurity and a new degenerate ground state appears with 
an intermediate coupling The FP corresponding to the ground state is characterized 

by non-zero residual entropy |^ and universal power low behavior in the impurity resistivity 
and the energy dependence of the scattering rate (-Rimp ~ T-*^/^, 1/r ~ tu^/^) [|]. Several 
normal state properties of Ce or y-based heavy fermion compounds like CeSn^ and YhCuAl 
can also be explained surprisingly well in terms of the 2-channel Kondo effect P]. 

Other extensively studied possible realizations of the two-channel Kondo model are pro- 
vided by fast two-level systems (TLS's) in metals These TLS's are formed by some 
heavy particles (HP's) tunneling between two neighboring sites and interacting strongly 
with the conduction electrons. Recent experimental and theoretical P,p!0| investigations 
confirm the conjecture that at low temperatures these TLS's can be properly described by 
an effective two-channel Kondo model, where the spatial motion of the HP is coupled to 
the angular momentum of the conduction electrons via an effective exchange coupling and 
the two degenerate channels correspond to the two real spin directions of the conduction 
electrons 



The TLS model sketched above is appropriate to describe tunneling centers in amor- 



phous metals |Tl|], however, it breaks down in systems like the narrow gap semiconductor 
Phi^xGcxTe or insulating Ki_xLixCl alloys |ll2| , p!3[| where the tunneling centers are formed 
by some substitutional impurities, and the HP is tunneling between 3, 6 or 8 equivalent posi- 
tions. Therefore the question arizes in a natural way, what is the low temperature behavior 
of an M-level system (MLS) which strongly interacts with the conduction electrons. 

To answer this question we consider a model previously introduced to describe the effect 



of the excited states in a TLS In the following we assume that the temperature (or 

the relevant energy scale) is low enough and thus the motion of the HP is restricted to the 
lowest lying M states corresponding to the M spatial positions of the HP: 
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M 

H,, = y: , (1) 

where hf creates a HP at site i and A*-^ is the tunnehng amphtude between positions i and 
j. We assume that no external stress is present and that the M positions are completely 
equivalent thus A** = 0. 

The most general two-particle interaction generated by the screened Coulomb interaction 
(or a pseudopotential) between the HP and the conduction electrons takes the form: 

^int= E KV^mh^nfO^^'mf . (2) 
i.j,n.m 

where creates a conduction electron with energy e, orbital quantum number n, and 
'flavor' index / (/ = 1,..,A^/). The orbital quantum number n can be thought of as the 
angular momenta of the electrons while for a real MLS the quantum numbers / = 1, 2 
correspond to the spin up and spin down conduction electrons and Nf = 2. For the sake 
of simplicity we also assume a constant density of states between the high and low energy 
cutoffs D and —D, qq, for all flavor numbers. Naturally, both the couplings and the 
tunneling amplitudes A*-^ are connected by the symmetry properties of the MLS which will 
be exploited later on. 



Following similar lines as in Ref. [T^ one can introduce the HP propagator Q^^{uj) and the 
HP-conduction electron vertex function F*-' (a;) in the usual way. However, calculating these 
functions in a perturbative way it turns out that they do not satisfy the usual multiplicative 
renormalization group (RG) equations. Therefore we consider the following generalized RG 
transformation (T = 0): 

g{u, v', A', D') = A g{u, V, A, D) A+ , 

r{uj, v', A', D') = [A+]-^r{u, V, A, D) A~^ , (3) 

where the matrix notations QqV^^ f , F^„ F, and A*-' A have been introduced, D' 
stands for the scaled bandwidth and A is an M x M matrix acting in the HP indices. Note 
that A = A{v' , A' , D' / D) is independent from the dynamical variable u. While for finite 



D/D' the matrix A has a rather comphcated structure for an infinitesimal change of D it 
can be chosen to be Hermitian and Eq. can be cast in the form of a scahng equation for 
the dimensionless couphngs v'^K 

In the following we assume that the relevant energy variable is cu, i.e., uj ^ \^'^^\,T. In 
this case the inverse HP propagator and the vertex functions can be expressed in the next 
to leading logarithmic order as 

id 

g^r^ = v'' - In ^ {[v'\ v''^] - Nfti{v"'v^^}v''^) , (4) 

where [ , ] denotes the commutator, the trace operator tr{...} is acting in the electronic 
indices, and a summation must be carried out over repeated indices. Plugging (|) into 
Eq. @ one can easily generate the scaling equations in a selfconsistent way: 



dx 
dv'^ 



— Nf 



where the scaling variable x = In (Dq/D) has been introduced. Since the scaling equations 
are rather complicated apart from some special cases they can be solved only numerically. 
However, to exploit the symmetry properties of the MLS it is useful first to introduce some 
site representation in the orbital indices of the conduction electrons. This can be achieved 
most simply by taking some linear combinations of the most strongly scattered angular 
momentum channels and hybridize them by using group theoretical methods. Working only 
with electron states directed to the impurity positions the f^-^'s become M^- dimensional 
tensors. However, the number of independent couplings is largely reduced by symmetry. 

A typical scaling of the norm of the dimensionless couplings, u = I] is shown in 

Fig. |I] (dashed line). The initial couplings have been estimated by using similar methods as 
in Refs. p|JT^ . Similarly to the multichannel Kondo problem both the infinite and the weak 
coupling FP's are unstable and the system scales to an intermediate strong coupling FP 
imj^,^. The Kondo energy can be identified as the crossover energy from the weak to strong 



coupling regimes: Tk = Dq e'^", Xc being the crossover value of the scaling parameter. For 
realistic initial parameters we find that this Kondo temperature can easily be found in the 
experimentally observable region, Tk ~ 1 — 10 i^'. 

In order to determine the properties of the MLS in the regime T, u; <C Tk one has to 
identify the FP the MLS scales to |T^. In the following we shall show that a MLS scales 
towards a FP, which - up to some potential scattering part - has the same structure as 



the SU{M) X SU{Nf) Coqblin-Scrieffer model ^7^. To prove this we first remark that the 
operators 5*-' J2k '^nm ^'^^ ^mn v^pp are invariant under scaling. Therefore v"^^ can be written 
as = + M^„, where the matrix M is a constant of motion depending on the initial 
parameters and — '^nn — 0- Then one can easily show that the right-hand side of 

Eq. (El) disappears provided 



1 fO'^ 0' 

— 



(6) 
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where the O^-^ 's are unitary equivalent to the generators of the SU{M) Lie algebra, 

with 01^^ ~ ^m^n ~ j^^^'^^nm- Similarly to the TLS problem beside the one in Eq. 
Eq. has an infinite number of FP's associated with different reducible and irreducible 
representations of the SU{M) Lie-algebra (J^). Our numerical simulations show, how- 
ever, that all the FP's which correspond to a representation different from the defin- 
ing one (0) are unstable. In Fig. |l| we show the scaling of the 'algebra coefficient' 
a = J2i,j,k,i \ \Nj[v^^ ,v^''] — Nf6^''v^^ + Nf6^^v^''\\, which measures in a natural way how well 
the fixed point (H) is approached by the w's. For T < Tk (i.e. for x > Xc) the coefficient a 
scales to zero. Thus below the Kondo temperature an MLS scales to the SU{M) x SU (Nf) 
Coqblin-Schrieffer model, which is a non- Fermi-liquid model and has a different scaling be- 
havior than the 2-channel Kondo model P,p!9|. (For a 6-state system, e.g., we expect an 
u^^^ behavior of the scattering rate.) 

To show that FP (^ is stable and to analyze its operator content we write the deviations 
of the couplings from their FP value in a form 



\ 



V 



(8) 



and linearize the scaling equations with respect to 6v'^K The couplings g^^ , f^, and /i*-' are 
M X M, M X oo, and oo x cxd matrices, respectively. Like the TLS case the linearized 
equations for g^^ , f^, and /i*-' decouple completely, 



dx N 



1- - M^.^ 



- 2n/ - 0*^tr{^^'^0^' + O^'^g''^} - tr{^*^0^'= + O^^'^^'^jO'^'} 



il±kk 



(9) 

(10) 
(11) 



and they can be solved exactly due to the simple structure of the O^-^'s. These linearized 
equations have an infinite number of zero modes; a finite number of them correspond to po- 
tential scattering while the others can be identified with the generators of the unitary trans- 
formations connecting the different possible M-dimensional subspaces where the SU{M) 
Lie-algebra is realized. These 0-modes can be shown, of course, to leave the Lie-algebra (|^ 
unaffected. All the other modes can be shown to be irrelevant, thus the FP (|^) is stable. 

The low energy properties of the MLS are determined by the operator content of the FP 
which is much richer and quite different from that of the simple Coqblin-Schrieffer model. 
A thorough analysis of Eq. (|lID shows that for M > 3 the leading irrelevant operators can 
be written as 



Oi 



V 



(12) 



C'^ 
{C^Y 

where the C^-^'s satisfy J2i{C^^ — CJI!^) = 0. These operators scale like ~ T'^' with A; = 
{M — 1)/Nf, and they describe scattering between channels which are not taken into account 
in the usual Coqblin-Schrieffer model. While they dominate the thermodynamical quantities 



like the specific heat, e.g., which scales as o 



imp 



T ' they do not contribute to the resistivity. 
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which scales hke ~ T^^' with \si = M/Nf, and is determined by subleading operators of the 
form 



O 



si 



r\-/ 



^ Q'i ^ 
y S'^ J 



(13) 



where the matrices Q'^ and S'^ satisfy E<5'' = = and = Q'^^. Note that the 

operators ([I^ ) do not exist in the TLS case (M = 2), and therefore the low-energy properties 



of a TLS can be completely described by the two-channel Kondo model |I01 . 

At this point we have to remark that since the FP couplings ■u*-' scale like the 
obtained scaling exponents can be considered as the first order estimates in a 1/Nf expansion, 
and they become exact in the Nj oo limit [p^OlpTSf (A^^ = M/Nf is e.g. the 1/A^j-expanded 
version of the exact exponent A^; = M/{Nf + M) for Nf > M l@,|19l). However, the 
expansion breaks down for Nf < M, and while for a physical MLS with Nf = 2 we expect 
similarly to other models that the fixed point structure remains the same, it remains 



an open question whether the new leading irrelevant operators (|T2|) survive in that physical 
limit or not. 

In order to decide whether the non-Fermi-liquid properties can be recovered in reality 
or not it is important to study the scaling of the splitting A*-' . As soon as the temperature 
(frequency) becomes smaller than the renormalized splitting |A*-^ (T, ci;)| the dynamics of the 
MLS is frozen out and the Kondo effect described above is stopped. Therefore in order to 
observe a non-Fermi-liquid behavior characteristic to the FP (^ we need \A'^^{Tx)\ ^ Tx- 
The scaling of the splitting parameter A^^ is shown in Fig. ^ (for the definition of A^^ see 
the inset in Fig. As one can see, for realistic model parameters A^^(Tft:)/A^^(Do) is as 
small as ~ 10^^, and therefore, even for very large splittings A^^ ~ lOOK the splitting is 
strongly reduced and since A^'^{Tk) ^ the MLS can get easily in the vicinity of the 
2-channel Coqblin-Schrieffer FP. 



As recently pointed out by Moustakas and Fisher for TLS's ||2^ multi-electron scattering 
becomes also relevant in the neighborhood of the two-channel Kondo fixed point. This 
process is also relevant in the present model in the vicinity of the NFL fixed point, however. 



similarly to the TLS case |2T[] it has a small initial amplitude, and being irrelevant in the 
weak coupling region this amplitude is even more significantly reduced during the first part 
of the scaling, Tk < D < Dq. Therefore, in the neighborhood of the FP it is always the 
splitting A*-' discussed above which provides the dominant mechanism to drive the system 
away the NFL FP and the multi-electron processes play a less important role. 

Concerning the experimental realizations there exist already some experiments on non- 
commutative 8-state systems in Pbi^xGcxTe alloys where the tunneling centers are formed 
by the relatively small Ge~^~^ ions [l^. While in this material a Kondo crossover has indeed 
been observed no non-Fermi-liquid behavior has been detected. However, in the experiments 
the inter-impurity interaction was quite strong which leads to the non-vanishing of the 
diagonal part A** of the tunneling matrix. Since, similarly to TLS's, the differences A** — 
A^^ are renormalized much less then the offdiagonal matrix elements [0, they can lead to 
the freezing out of some impurity states corresponding to transitions from the SU{M) to 
SU{M') models with M' < M, finally, most probably, reaching the trivial M' = 1 model. 
Furthermore, in this material there is a strong spin-orbit scattering, which leads to a strong 
cross scattering between the spin up and spin down electron channels and drives the system 
to the Nf = 1 Fermi-liquid FP instead of the Nf = 2 non-Fermi-liquid model. Thus, for this 
material even if the inter-impurity interactions were small we would expect a Fermi liquid 
behavior. Therefore it would be very interesting from the experimental side to find some 
similar alloys with small spin- orbit interaction, where the spin symmetry would guarantee 
the non-Fermi-liquid properties. To observe the non-Fermi liquid FP one could apply a 
pressure on the sample, which is an appropriate tool for tuning into the region 1-lOK 
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FIGURES 

FIG. 1. Scaling of the norm of the dimensionless couphngs, u = J2 Ik*'' 1 1 (dashed hne), and of 
the algebra coefficient a (continuous line) for a 6-state system with Nf = 2. The initial couplings 
have been chosen to be vH = 0.8, v^l = 0.2, = 0.1, vjl = 0.05, = 0.03, = v% = 0.0005, 
vlf = 0.005, vlf = 0.003, Vqi = vj^ = 0.0005. The other non-zero couplings have been generated 
by symmetry transformations. 

FIG. 2. Scaling of the dimensionless hopping amplitude, A^'^/Dq for the same 6-state system 
as in Fig. 1. Inset: Numbering of the sites of the 6-state system. 
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